Accepted for publication in The Astrophysical Journal 

Preprint typeset using KTgX style emulateapj v. 1 1/10/09 



THE DETECTABILITY OF TRANSIT DEPTH VARIATIONS DUE TO 
EXOPLANETARY OBLATENESS AND SPIN PRECESSION 

Joshua A. Carter & Joshua N. Winn 

Kavli Institute for Astrophysics and Space Research, 

Massachusetts Institute of Technology, Cambridge, MA 02139 

carterja@mit.edu 

Accepted for publication in The Astrophysical Journal 



o 

X. 



6 



> 

o 
o 



% 



ABSTRACT 

Knowledge of an exoplanet's oblateness and obliquity would give clues about its formation and internal 
structure. In principle, a light curve of a transiting planet bears information about the planet's shape, but 
previous work has shown that the oblateness-induced signal will be extremely difficult to detect. Here we 
investigate the potentially larger signals due to planetary spin precession. The most readily detectable effects 
are transit depth variations (T(5V) in a sequence of light curves. For a planet as oblate as Jupiter or Saturn, 
the transit depth will undergo fractional variations of order 1%. The most promising systems are those with 
orbital periods of approximately 15-30 days, which is short enough for the precession period to be less than 
about 40 years, and long enough to avoid spin-down due to tidal friction. The detectability of the T^V signal 
would be enhanced by moons (which would decrease the precession period) or planetary rings (which would 
increase the amplitude). The Kepler mission should find several planets for which precession-induced TSW 
signals will be detectable. Due to modeling degeneracies, Kepler photometry would yield only a lower bound 
on oblateness. The degeneracy could be lifted by observing the oblateness-induced asymmetry in at least one 
transit light curve, or by making assumptions about the planetary interior. 
Subject headings: stars: planetary systems — techniques: photometric 



1. INTRODUCTION 

Measuring the oblateness of exoplanets would further our 
understanding of planetary formation, rotation, and internal 
structure. One possible measurement technique relies on 
the differences between the transit light curve of a spheri- 
cal planet and an oblate planet with the same sky-projected 
area (Seager & Hui 2002, Barnes & Fortney 2003). However, 
the differences are minuscule, of order 200 parts per million 
(ppm) for a planet as oblate as Saturn, and 2 ppm for a "hot 
Jupiter" whose spin rate has been slowed by tidal friction into 
synchronization with its orbit. In a previous paper we showed 
that with this technique, even the best available light curves 
are only barely sufficient to rule out a Saturn-like oblateness 
(Carter & Winn 2009). Those results pertained to the planet 
HD 189733b, for which one would expect spin-orbit synchro- 
nization, and consequently the theoretical oblateness was an 
order of magnitude below the empirical upper limit. 

In that work we also pointed out the potentially observable 
effects of a phenomenon that had been previously overlooked: 
the precession of the planet's rotation axis. Precession of 
an oblate planet causes the sky-projected area of the planet 
to change over time, thereby causing gradual changes in the 
depth and duration of transits. In this paper we investigate 
the observable manifestations of spin precession in a broader 
context. Section 2 describes the characteristics of the signal 
in terms of the properties of the star, planet, and orbit. Section 
3 discusses the range of orbital periods for which the signal 
is most readily detectable, considering the timescales for 
precession and spin-orbit synchronization. Section 4 presents 
simulated results for a specific case, a Saturn-like planet 
observed by the Kepler satellite. Section 5 summarizes and 
discusses the results. 



2. CHARACTERISTICS OF THE SIGNAL 

We model the planet as an oblate spheroid, illustrated in 
Figure ([1]). The oblateness (or flatness) parameter is defined 

as 
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where /?eq and Rpoi are the equatorial and polar radii, respec- 
tively. For rotationally-induced oblateness, a good approxi- 
mation is 
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where Mp is the mass of the planet and J2 is the planet's 
zonal quadrupole moment (Murray & Dermott 2000, Hubbard 
1984). This approximation is valid for all the Solar System 
planets for which /, J2 and P,ot have been measured precisely. 
The obliquity is the angle 9 between the polar axis and the 
orbital axis. The angle (f> specifies the direction of the projec- 
tion of the polar axis onto the orbital plane. For a uniformly 
precessing planet, the case considered here, is a linear func- 
tion of time. 
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Accurate calculations of the transit light curve of an oblate 
spheroid, including the effects of stellar limb darkening, have 
been presented by Seager & Hui (2002), Barnes & Fort- 
ney (2003), and Carter & Winn (2009). In this paper we are 
not interested in high accuracy or in the slight differences be- 
tween the light curve an oblate planet and a spherical planet. 
Instead we are interested in the order of magnitude of the vari- 
ations in the transit depth and duration due to the changes in 
the precessing planet's sky-projected figure. 




Fig. 1 . — Geometry of transits by an oblate spheroid. The planet's orientation is specified by the spin-orbit obliquity d = cos"' s ■ n, the angle cp between the 
projection of the spin axis onto the orbital plane and the vector connecting the centers of the star and planet at midtransit, and the orbital inclination i relative to 
the sky plane. 



In the absence of limb darkening, the depth 5 is approxi- 
mately the areal ratio between the sky projection of the oblate 
spheroid and the stellar disk, 
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where k = Rf^^/Rj, is the planet-to-star radius ratio, / is the 
orbital inclination with respect to the sky plane and e is the 
ellipticity. 
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A derivation of this expression is given in the Appendix. Fig- 
ure|2]shows the fractional amplitude of the transit depth varia- 
tions for the case i = 90°, as a function of/ and 9. For Saturn- 
like values of oblateness and obliquity, the depth variations 
would be a few percent. 

The transit duration will also vary over the precession pe- 
riod, due to the changing dimension R\\ of the planet's sky 
projection in the direction of orbital motion: 
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where 
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are the quantities describing the oblateness and obliquity of 
the projected exoplanet's shape. Derivations of these expres- 
sions are also given in the Appendix. For a circular orbit, 
the ingress/egress duration (first to second contact, or third to 
fourth contact) is approximately 
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where forb is the orbital period, a is the orbital distance, and 
b = a cos i/R^, is the normalized impact parameter. The full 
transit duration (first to fourth contact) is approximately (e.g.. 



Seager & Mallen-Ornelas 2003) 
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These approximations are valid as long as the transit is not too 
close to grazing. The fractional amplitude of the r variations 
(TtV) is comparable to that of depth variations (T(5V). The 
amplitude of transit full-duration variations (TDV) depends 
on k and b, in addition to / and 6, and therefore cannot be 
summarized in a single contour plot such as Figure |2] 

3. THE OPTIMAL ORBITAL DISTANCE 

In this section we suppose that the planet's precession is 
caused exclusively by the gravitational torque from the star. In 
order for the detection of precession-induced T(5Vs or TDVs 
to be feasible, the planet must be close enough to the star for 
precession to produce observable effects in a human lifetime. 
However, if the planet is too close to the star, then tidal dissi- 
pation should slow down the planet's rotation until it is syn- 
chronized with the orbital period, and drive the obliquity to 
zero, which would cause the signal to be undetectable. Hence, 
we must ask if there is a range of distances from the star that 
is close enough for rapid precession, and yet far enough to 
avoid spin-orbit synchronization. 

The spin precession period for a planet on a fixed circular 
orbit is given by 
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(Ward 1975), where Piot is the planet's rotation period and C is 
its moment of inertia divided by MpRl . The numerical scal- 
ing of 13.5 for C/J2 is the estimated value for Saturn (Ward & 
Hamilton 2004). According to this expression, orbital periods 
shorter than Poib ^30 days will lead to rapid enough preces- 
sion to be observed over decadal timescales, depending on 
the planet's obliquity and internal structure. In Fig. [3] the 
thick solid line shows the spin precession period as a function 
of orbital distance, for an exoplanet with the same Pjot, C, J2, 
and 9 as Saturn. (The thin solid lines show the more rapid pre- 
cession rates produced by hypothetical planetary satellites, as 
discussed in §|5]) 
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Fig. 2. — The fractional amplitude of transit depth variations (T5Vs), shown as a function of oblateness (/) and obliquity (0) assuming the orbit is perpendicular 
to the sky plane ((' = 90°). For reference, the astrological symbols show the shape parameters of Solar System planets: ^ - Mercury, 9 - Venus, ffi - Earth, d" - 
Mars, % - Jupiter, f} — Saturn, 5 - Uranus, ^ - Neptune. 
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Fig. 3. — The optimal orbital distance. The solid curves refer to the axis 
on the left; they show the calculated spin precession period as a function of 
orbital distance for a Saturn-analog planet, with no satellites (thick line) and 
with a Titan-like satellite (thin lines) orbiting at a specified fraction of a Hill 
radius. The dashed curve refers to the axis on the right; it shows the calculated 
timescale for tidal spin-orbit synchronization. The vertical dotted line marks 
the orbital distance for which the spin-down time is 1 Gyr 

The approximate timescale for tidal spin-orbit synchroniza- 
3n is 
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(Goldreich & Soter 1966), where Qp is the planet's tidal dis- 
sipation factor and Pyotj is its initial rotation period, both of 
which are highly uncertain. 

Based upon Solar System constraints, Qp is thought to be 
in the range 10-500 for terrestrial planets, and 10'*^"^^ for 
gas or ice giant planets (Goldreich & Soter 1966, Peale et 
al. 1980, Yoder 1995, MardHng & Lin 2004, Ogilvie & Lin 
2004, Jackson et al. 2008), although the results are strongly 



model-dependent and also dependent on the frequency of tidal 
oscillations. In setting the scale parameters in Eqn. (fT2l i we 
adopted Qp = 10^^, on the high end of current estimates, giv- 
ing the longest (most favorable) synchronization timescale. 

The primordial spin period might be expected to be near the 
rotational breakup limit, which is w3 hr for Jupiter (Murray & 
Dermott 2000), although the effects of planetary contraction 
and disk-planet interactions should also be considered. For 
Eqn. (fT2l i we used a current rotation period of 10 hr (similar 
to Jupiter and Saturn) and an initial period of 5 hr. 

For ^s ^ 1 Gyr, it is reasonable to hope that the planet 
has not yet been tidally spun down. Thus, with reference to 
Eqns. (fllT i and ( fT2] i. the "sweet spot" for observing the effects 
of spin precession on the transit parameters is at Porb ^ 15 d, 
which is long enough to allow for rapid rotation, and short 
enough to allow for rapid precession. 

Of course the identification of a single optimal period is a 
simplification. The existence and observability of T^Vs and 
TDVs depends on the particular mass and radius of the planet 
under consideration and the observed age of the star, as well 
as the parameters relating to the planet's tidal dissipation, in- 
ternal constitution, and initial spin period. The dashed line in 
Figure [3] shows the dependence of Tj on orbital distance for 
the particular case of a Saturn analog. Table [T] gives some 
numerical results for precession periods and synchronization 
timescales for hypothetical close-in planets with properties 
similar to Solar System planets. If we require Ppiec < 40 yr 
and Ts > 1 Gyr for observability, then from Table [T] we see 
that the signal is potentially observable for Jupiter and Saturn 
analogs. For analogs of Uranus and Neptune, we expect the 
signal to be unobservable unless the spin precession is made 
more rapid by the presence of large satellites (see § |5]l. For 
Earthlike planets the signal seems unlikely to be observable 
because of strong tidal dissipation. 

4. A SPECIFIC EXAMPLE 

In this section we examine the particular case of a hypo- 
thetical, ringless, moonless, Saturn-like transiting exoplanet 
(/ = 0.1, 6 = 27°) in a circular orbit around a Sun-like star 



TABLE 1 
Relevant timescales for the observability of planetary spin precession 
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C//2 
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r, [Gyr] 


Observable Range in Pc-i, [d] 


Jupiter 


5.0 


10.0 


10^5 


0.26 


15.4 


13.7 


14.3 


24.3 


8.5 


14.3 - 24.3 


Saturn 


5.0 


11.0 


106.5 


0.22 


13.5 


15.7 


17.1 


27.3 


6.5 


17.1-27.3 


Uranus 


8.0 


17.0 


106.5 


0.22 


67.0 


45.3 


16.3 


15.3 


0.78 


None" 


Neptune 


8.0 


16.0 


106.5 


0.23 


68.0 


43.8 


15.3 


14.7 


0.83 


None" 


Earth 


12.0 


24.0 


102.0 


0.34 


314 


13,100 


151 


8.40 
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None" 



References. — Murray & Dermott (2000), Hubbard (1984), Yoder (1995), Ward & Hamilton (2004) 
' Assuming no satellites 



with Porb = 17 d. This orbit is near the "sweet spot," giving 
^prec = 17 yr and Ts = 1 Gyr. 

4.1. Likelihood of discovery 

First we must ask how likely it is that such a planet will be 
discovered. Giant planets with periods between 15-30 d are 
already known to exist from Doppler surveys. At the time of 
writing, the exoplanets . org database has 8 such planets. 
Indeed one of them is already known to transit (HD 17156b; 
Fischer et al. 2007), although in that case tidal effects may 
have slowed the planet's rotation because the orbit is highly 
eccentric and the pericenter distance is small. None of the 
other 7 planets is known to transit, and the probability that at 
least one of them transits is approximately 23%. Given these 
facts and the recent acceleration in the discovery rate using 
the Doppler method, it would seem likely that a transiting gas 
giant with Porb = 15-30 d will be discovered in the near future. 

Of particular interest are the prospects for the Kepler satel- 
lite mission (Borucki et al. 2010). Kepler searches for tran- 
siting planets by keeping wlO^ stars under nearly continuous 
photometric surveillance for at least 3.5 yr and possibly as 
long as 6 yr. To estimate the fraction of the target stars that has 
a suitable planet, we used the power-law formulas given by 
Gumming et al. (2008) for the abundance of planets of a given 
mass and period, which were derived from data from the Keck 
Planet Search. The integrated value of (abundance x transit 
probability) over the range Porb =15-30 d, Mp = 0.2-2 Mjup 
is approximately 0.019%. Therefore, among the 5 x lO'* Sun- 
like stars in the Kepler field, we expect wlO transiting giant 
planets with PqA = 15-30 d. 

4.2. Expected signal 

Figure |4] shows the expected time variations in S, t, and 
Tfuii for our hypothetical close-in Saturn. Compared to the 
TSV signal, the TDV signal has a smaller amplitude, and de- 
pends on a larger number of parameters. The TrV signal has 
a similar amplitude to the TSV signal, but t cannot be mea- 
sured as precisely as S because of the relatively short duration 
of the partial transit phases. Therefore we expect the T5Y 
signal to provide the best constraints on the planetary shape 
parameters. 

Depth variations of a few percent should be detectable. 
Ground-based observations of individual transits have already 
allowed S to be measured to within 1% [see, e.g., Gillon et 
al. (2009), Johnson et al. (2009), Winn et al. (2009)]. To es- 
timate the precision with which Kepler can measure transit 
depth variations, we simulated a transit light curve and aver- 
aged it into 30 min bins, corresponding to Kepler's standard 
time sampling. Then we added white noise with a standard 
deviation of 95 ppm, as appropriate for a typical Kepler tar- 



get (KIC r mag = 13), based on the early results of Welsh et 
al. (2010) and Latham et al. (2010). By fitting a parameter- 
ized model to the light curve we found that the transit depth 
was recovered to within 0.6%. Over a 6 yr mission, Kepler 
would observe «130 transits with this precision. 

4.3. Parameters and degeneracies 

Next we discuss the "inverse problem" of inferring the plan- 
etary shape and precession parameters from the observed TJV 
signal. For / = 90°, the TJV signal is periodic with period 
PpKc/2. For / ^^'90° the TSV period is actually Pprec, but the 
odd and even maxima have only slightly different amplitudes. 
Hence the period of the TSY signal will reveal Pprec- 

Three other characteristics of the TSY signal will help to 
pin down the model parameters. First, the minimum observed 
transit depth will reveal the parameter combination 
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Second, the full range of transit depths determines a combi- 
nation of / and 6, 
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where A/ = 7r/2 - / is always small and can be measured pre- 
cisely through transit photometry. Third, the phases at which 
the maxima occur will determine 0o- 

For A/ y there are two other informative observables, at 
least in principle. First, the phases at which the minima occur 
depends slightly on A/ and d\ 
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heights between odd and even maxima. 
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but these effects are very small in practice. For our Saturn 
analog, assuming / = 88.97° (transit impact parameter = 0.5), 
the difference in the maxima is smaller than 0.25%. 

For A/ y there are a total of 6 observables listed above, 
and it would be possible in principle (with arbitrarily precise 
data) to determine all 6 model parameters k, f, 6, i, Pprec and 
(f)o. However, using our simulated Kepler Ught curve we find 
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Fig. 4. — Variations in tlie transit light curve due to an oblate, oblique, precessing exoplanet. Plotted are the transit depth (5), total duration (rf„ii) and ingress 
duration (r) fractional variations (T<5V , TDV, and TrV, respectively) that are expected for a uniformly precessing Saturn-like planet around a Sun-like star. The 
time scale is based on the assumption Port = 17.1 days. 



that / and 9 cannot be determined independently, although 
Eqn. (O could be used to place a lower bound on /. The 
degeneracy is illustrated in Fig. (|5]l. 

To break the parameter degeneracy, one possibility is to ar- 
range for high-cadence, high-precision observations of at least 
one transit, seeking the slight oblateness-induced anomalies 
that were described by Seager & Hui (2002) and Barnes & 
Fortney (2003). Observations of a single light curve would 
lead to constraints on the sky-projected oblateness and obliq- 
uity [see Eqns. d?) and (O], which, together with the TSY 
signal, would uniquely determine / and 9. The best time to 
schedule such observations would be near a minimum of the 
T(5V curve, when the light curve anomalies are largest. 

This would be a challenging task, as the amplitude of 
the differences between the actual light curve and the best- 
fitting model of the transit of a spherical planet would be 
<100 ppm. For this specific example, even Kepler photom- 
etry (with 1 min cadence) would be insufficient to detect the 
anomalies in a single transit. By observing 10 transits near 
the minimum of the TSW signal (during which time the sky- 
projected quantities are constant to within 1%), Kepler could 
detect the signatures of oblateness and obliquity at the Icr 
level, but the resulting constraints would be weaker than the 
constraints determined by an analysis of the T(5V signal. A 
significantly larger planet, or brighter host star, would be re- 
quired for meaningful constraints. 

Another possibility is to enforce additional physically- 
motivated relationships between parameters. In particular, we 
have already shown that the precession period is a function of 
k, /, 9, C and Ji [by combining Eqn. (|2]l with Eqn. (fTTTil. A 
further condition can be imposed on /, C and J2, such as the 
Darwin-Radau approximation for planets in hydrostatic equi- 
librium (Murray & Dermott 2000): 



(17) 



Following this path, there are 8 model parameters {k, f, 9, 
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among them. However there are only 5 quantities that are 
well determined from the photometric data [Pprec, (po, /sin 9, 
k^{l -/), /], leaving us still short by one observable or con- 
straint from being able to determine all the parameters. For 
example, if one were willing to assume C = 0.23, then we find 
from our simulated Kepler data that /, 9, J2 and Ppi-ec can be 
recovered with a precision of about 10%, and k is recovered 
within 1%. 

Assigning C a specific value is unrealistic, but for realis- 
tic planets one expects C to be smaller than 0.4 (Murray & 
Dermott 2000). We repeated the analysis of our hypotheti- 
cal T(5V signal, allowing C to be a free parameter restricted to 
that range, with a uniform priorOl In effect we averaged the 
results over a range of C deemed to be physically plausible. 
As might be expected, a strong degeneracy was observed be- 
tween / and 9, as seen in Fig. ©• However, we were still able 
to determine J2 to within 10%, and k to within 2%. 

5. DISCUSSION 

In this paper we have investigated the observability of 
changes to transit light curves resulting from the spin pre- 
cession of an oblate, oblique exoplanet. The most readily 
detectable signal is the T(5V signal, the variation in transit 
depth due to the changing area of the planetary silhouette. 
The planets that seem most likely to exhibit detectable effects 
are those with periods between 15-30 days (around Sun-like 
stars), which is short enough for precession periods to be 40 yr 
or less, and long enough to hope that tidal spin-orbit synchro- 
nization has not taken place. 

It is also important to consider other physical processes that 
could give rise to T(5V signals, and which might confound 
the interpretation of the data. Starspots and other types of 
stellar variability can produce transit depth variations. These 
can be recognized and taken into account by monitoring the 
star outside of transits, as is done automatically by the Kepler 

' We also required J2 > 0, which corresponds to C > 0.133 according to 
the Darwin-Radau relation (Eqn. llVt . 
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Fig. 5. — Hypothetical constraints on / and 6, as determined by analyzing 
a simulated T<5V signal measured by Kepler. The signal was based on transits 
of a Saturn analog on a 17.1-day orbit around a Sun-like star with KIC r= 13 
mag. The light gray region marks the 95% confidence region, using a model 
that allow the parameters {k,f,6,i,PpTcc,4'o} to be independent free parame- 
ters. In this case only the product fsir? 9 is well determined; the dashed line 
satisfies /sin^ 9 = constant. The dark gray region is the 95% confidence re- 
gion, using a model that enforces the expected physical relationship between 
T(5V amplitude and period, as well as the Darwin-Radau relation and a re- 
stricted range for the normalized moment of inertia C. See §|4]for details. 
The dotted lines mark the "true" values of / and 9 that were used to generate 
the simulated T(5V signal. 

satellite. We are not aware of any atmospheric phenomena 
associated with the exoplanet that would result in variations 
in the projected area at the 1 % level, mimicking those due to 
uniform spin precession. 

However, we can think of two plausible phenomena that 
would affect the T5N signal: moons and rings. If the planet 
has any moons then the precession period may be shorter than 
we have calculated. Satellites provide more leverage for the 
star to torque the exoplanetary system. Mathematically, satel- 
lites augment the effective values of J2 and C. Following 
Ward & Hamilton (2004) we may write the enhanced values 
as J2 + j and C + c, where 



c = 



V^ ?M; 



at 
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where m,, a,, /,■ and p'^^ are the satelUtes' masses, orbital radii, 
orbital inclinations (relative to the planetary equator), and or- 
bital periods. 

For example, in the Saturnian system j/Ji ~ 3.2 while 
c/C w 0.01 such that (C + c)/(72 + ;) « 3.2 (as compared to 
C/72 = 13.5). Titan alone is responsible for about 90% of j 
and c, shortening Saturn's precession period by a factor of 
four relative to a satellite-free SaturnQ Fig. (O shows the ef- 
fect of Titan analogs at various distances around the Saturn 
analog considered in this paper. The distances are expressed 



as fractions of the Hill radius rn- This effect might be used 
to implicate the presence of exomoons, if a T^V signal were 
observed and found to correspond to an effective J2 too large 
to be plausibly attributed to the planet alone. 
Ring systems would increase the amplitude of the signal. 



^ For Titan: a,/Req ■ 
& Dermott (2000). 



: 20.2, mi/Mp = 2.3 X lO-*, pjj^^/Proc = 38.3; Muiray 



while leaving the period unchanged. Optically thick rings 
that lie within the equatorial plane of the planet would in- 
crease the fractional variation of the sky-projected area, as 
the planet precesses. In contrast to exomoons, ring systems 
have little mass, and would not significantly reduce the pre- 
cessional period. This is the situation in the Saturnian sys- 
tem (Ward & Hamilton 2004). It is not certain that rings 
could exist around planets having orbital periods between 15 
and 30 days. Rings around planets with orbital periods less 
than Poib ~ 15 days would likely be short-lived as a result 
of Poynting-Roberston drag amongst other destructive effects 
(Barnes & Fortney 2004). Also, we would be unlikely to find 
rings composed of water ice (e.g. Saturn's rings) around plan- 
ets whose orbits are interior to the snow line (< 1 AU). Nev- 
ertheless, rings of other compositions may exist. 

In short, rings and moons would each affect the observed 
T(5V signal, and in complementary ways. This may introduce 
some ambiguity in the estimation of the shape parameters of 
the rotating planet, but may also allow the rings and moons to 
be detectable, issues that we leave for future work. 

For simplicity we have considered only circular orbits. 
Planets on eccentric orbits will undergo apsidal precession 
and nodal precession, which will result in time variable stel- 
lar torques on the planet and consequent modifications to 
the spin-axis precession. For Saturn or Jupiter analogs at 
Porh ^15 days, the apsidal and nodal precession periods are 
>10^ yr and therefore likely to be irrelevant (Ragozzine & 
Wolf 2009). 

In addition, we restricted our attention to the simplest case 
of uniform spin precession, but in reality the perturbations 
from other bodies may cause the spin axis to perform a more 
complex ballet. For example, Mars's spin axis tumbles chaoti- 
cally (Touma & Wisdom 1992). Saturn's moon Titan causes a 
700 yr modulation of Saturn's spin precession frequency, due 
to its inclined orbit. Furthermore, Saturn's spin axis may be 
trapped in a resonance with Neptune's orbit, causing it to li- 
brate with an angular amplitude of >31° as it circulates about 
the second Cassini state (Ward & Hamilton 2004). These 
effects would be manifested as additional time dependences 
("noise") in the T(5V signal. The effects are impossible to fore- 
cast for exoplanets, depending as they do on the existence of 
other bodies and any resonances that may occur. 

We thank Dan Fabrycky, Darin Ragozzine, and members 
of the MIT exoplanet discussion group, for helpful conversa- 
tions. We also thank an anonymous referee for helpful com- 
ments on an earlier draft of this manuscript. 



APPENDIX 
THE PROJECTION OF A SPHEROIDAL PLANET ONTO THE PLANE OF THE SKY 

We model the planet as an oblate spheroid, illustrated in Figure ([T]). The sky plane projected shadow of the spheroidal planet is 
bounded by an ellipse. The polar axis, s, of the spheroid is tilted by the obliquity angle 6* > from the orbital axis h and by the 
angle 9' > Q from the axis h' which is perpendicular to the plane that is perpendicular to the sky plane. Let y be the axis which 
lies along the line connecting the center of mass and the observer Let i be the axis in the orbital plane whose projection into the 



plane perpendicular to the sky plane lies along y. The angle (f) is the angle through which s is rotated about n and 0' is the angle 
through which s is rotated about n' such that = 0' = corresponds to the orientation in which the spheroid is tipped towards the 
observer (such that s is coplanar with both i and y). It follows that 

(Al) 

(A2) 
(A3) 
(A4) 
i-3' = n cos/-i (l-sin/) (A5) 

where / = sin"' x-y is the inclination of the orbital plane to the sky plane and where we have assumed the distance from planet 
to star is much less than the distance from planet to observer. The following relationship between unprimed angles (measured 
relative to the orbital plane) and primed angles (measured relative to the sky plane) may be derived from the above equations: 

cos0'sin6'' = sin/ cos0sin6' + cos/ cos6'. ( A6) 

The polar axis of the spheroid is inclined relative to the plane of the sky by the angle 9" where 



cos u = s 


n 


cos9' = s 


h 


cos<j)sin9 = s 


X 


cos(f>'sin9' = s 


y 



cos 



2 nil 



= l-cosVsin^6l' (A7) 

= 1 - (sin / cos (/) sin + cos / cos 9) ( A8) 

If p and q are the orthonormal axes in the sky plane that are aligned with the axes of the ellipse bounding the planet's shadow, 
then a point {y,p,q) on the surface of the spheroid satisfies 

/ ^' \2 

(A9) 

where 



rv'i2 1 ,,2 1 M^ 


= < 


(y)^p^[,_f) 


y' =ycos9"-qsm 


9" 


q' = ysin9" + qcos 


9", 



7?eq is the equatorial radius and / is the oblateness parameter. 

The analytic description of the closed curve bounding the projection of the planet on the sky plane may be found by solving 
for the non-degenerate solutions of yip, q) in Eqn. (IA9b . Non-degenerate solutions correspond to points in the sky plane where 
rays parallel with y intersect at exactly one point on the planetary surface and thus define the boundary of the shadow as seen by 
the observer (who is located at a distant point along y). As advertised, it may be shown that the collection of these points satisfy 
the equation of an ellipse, 

sm 6'" + (l-/)2cos2 6'" ^ 
having major and minor axis lengths, A and B, satisfying 

A=/?eq (All) 



B = R^^^ sin^ 6"' + ( 1 - fY cos2 9" 



= /?eqVl-e2cOs2 6'" 



= ^eqV 1" 



l-(sin/ COS (/) sin 6* + cos/ cos 6*) 



(A12) 



where e = ^J\-(\-fY. 

The minor axis of the ellipse lies along the projection of the polar axis onto the sky plane. The angle between the major axis of 
the ellipse and the direction of orbital motion is given as 

_, sinfiitan^ ,. ,^^ 

9i_ = tan ' -—^ —^. (A13) 

cos tan 9 cos ; - sm ; 

and the oblateness parameter of the ellipse, f± = (A- B)/A, is given as 

/j^ = 1- W l-e^i l-[sin0cos(/)sin/ + cos6'cos/]"l. (A14) 

The boundary of the ellipse may be defined relative to the ellipse center via the parameter equations 

X{s) = Req [cos.scos0^-(l-/x)sin.ssin0^] (A15) 

y(i) = /?eq [cos5sin0_L + (l-/^)sin5COS0j_] (A16) 
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for s e [0, 2tt] and where X is the coordinate along the direction of orbital motion. The radius, /?|| , of the projection of the ellipse 
onto the direction of orbital motion is equal to the maximum ofX(s). This maximum occurs for s^ax satisfying 



tan.Sn,ax = -(l-/L)tant 



for which 



The areal ratio of the spheroidal planet's sky projection and the stellar disk is 

, ttAB 

= 



(A17) 



(A18) 



ttRI 



l-e2 



l-(sin/ cos (/) sin 6' + cos; cos 6*) 



(A19) 



where k = Req/R^, 
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